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Abstract 

In the standard cosmological framework of the Oth-order FLRW metric and the use of perfect 
fluids in the stress-energy tensor, dark energy with an equation-of-state parameter w < — 1 (known 
as phantom dark energy) implies negative kinetic energy and vacuum instability when modeled as 
a scalar field. However, the accepted values for present-day w from Planck and WMAP9 include 
a significant range of values less than —1. We find that it is not as obvious as one might think 
that phantom dark energy has negative kinetic energy categorically. Analogously, we find that 
field models of quintessence dark energy ( w^ > — 1) do not necessarily have positive kinetic energy 
categorically. Staying within the confines of observational constraints and general relativity, for 
which there is good experimental validation, we consider a few reasonable departures from the 
standard Oth-order framework in an attempt to see if negative kinetic energy can be avoided in 
these settings despite an apparent w < — 1. We consider a more accurate description of the universe 
through the perturbing of the isotropic and homogeneous FLRW metric and the components of 
the stress-energy tensor, and we consider dynamic w and primordial isocurvature and adiabatic 
perturbations. We find that phantom dark energy does not necessarily have negative kinetic energy 
for all relevant length scales at all times, and we also find that, by the same token, quintessence 
dark energy does not necessarily have positive kinetic energy for all relevant length scales at all 
times. 


* kludwick@virginia.edu 


1 




Introduction 


A recent milestone in observational cosmology happened when the High-z Supernova 
Search Team in 1998 [I] and the Supernova Cosmology Project in 1999 [2j published obser¬ 
vations of the emission spectra of Type la supernovae indicating that the universe’s rate of 
outward expansion is increasing. Galaxy surveys and the late-time integrated Sachs-Wolfe 
effect also give evidence for the universe’s acceleration. Thus, ’’dark energy” was proposed as 
the pervasive energy in the universe necessary to produce the outward force that causes this 
acceleration, which has been observationally tested and vetted since its discovery. The 2011 
Nobel Prize in Physics was awarded to Schmidt, Riess, and Perlmutter for their pioneering 
work leading to the discovery of dark energy. The present-day equation-of-state parameter 
w from the equation of state most frequently tested by cosmological probes, p = wp with 
constant w, assuming a flat universe and a perfect fluid representing dark energy, has been 
constrained by Planck in early 2015 to be w = —1.006 ± 0.045 [3], and Planck’s 2013 value 
was w = — 1.131q lo 0 • The value from the Nine-Year Wilkinson Microwave Anisotropy 
Probe (WMAP9), combining data from WMAP, the cosmic microwave background (CMB), 
baryonic acoustic oscillations (BAO), supernova measurements, and Ho measurements, is 
w = —1.084 ± 0.063 [5]. From these reported values, the prospect of w < —1 is clearly a 
distinct possibility, and under other assumptions (such as a spatially curved universe), the 
window reported for w does not always include the value for the cosmological constant (CC) 
model, w — — 1. 

However, dark energy modeled as a perfect fluid with w < — 1 leads to a field theory 
with negative kinetic energy (a ghost field theory), which implies vacuum instability. Either 
the phantom ghost has positive density and violates unitarity, rendering it unphysical, or 
unitarity is satisfied and the density is negative, which leads to vacuum instability [6]. This 
phantom dark energy with a wrong-sign kinetic term described as an effective field theory 
may be able to make this instability unobservable, but not without great difficulty and 
perhaps sacrifice of well-accepted physical principles 00- 

One deduces the ghost nature of phantom dark energy from w < — 1 within the standard 
cosmological framework of the Oth-order Friedmann-Lemaitre-Robertson-Walker (FLRW) 
metric with the use of perfect fluids in the stress-energy tensor, but the condition for negative 
kinetic energy is different for different frameworks. In this work, given that our universe 
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is not perfectly isotropic and homogeneous, we examine the possibility for positive kinetic 
energy with w < —1 in light of first-order perturbations to the FLRW metric and the 
components of the stress-energy tensor. We first consider constant w and then dynamic 
w. We also consider primordial isocurvature along with adiabatic perturbations. We then 
consider the possibility of negative kinetic energy for w > — 1 via the inclusion of cosmological 
perturbations. Although there are many alternative frameworks one may study (inclusion 
of spatial curvature, dynamic sound speed c s , vector field dark energy, models with coupled 
dark energy and dark matter, modified gravity, different metrics, dark energy as an imperfect 
fluid, quantum corrections), we focus on the aforementioned manifestations of dark energy in 
this paper and leave these others to future work. There is abundant experimental verification 
of general relativity, so we study the condition on the kinetic energy of phantom dark energy 
within the confines of general relativity and observational constraints from cosmological 
probes. 


Phantom Dark Energy 


Consider the Einstein-Hilbert action for general relativity with a complex scalar field 
(c=l): 


S = j d A x\f^g 




+ Smi 


( 1 ) 


where the first term is the usual contribution to the Einstein tensor, the second and third 
terms are the contribution to the scalar field dark energy, and S m is the action for the rest 
of the components of the stress-energy tensor Minimizing the action leads to Einstein’s 
equation, 

Rixu - l^Rg^u = 8vrG(T^[0] + T^lm]), (2) 

where T^ u [(p} = -2 ^ + g^C#. 

Assuming dark energy is spatially homogeneous as a perfect fluid, the density p$ and 
pressure P ( p for the scalar field are 


P<t> 


2 a 2 


+ VM), p* 


2 a 2 


-VM). 


( 3 ) 


We used the flat FLRW metric 


ds 2 — a 2 (r ) —dr 2 + dx l dxi , 


( 4 ) 
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and • represents differentiation with respect to r. 

The kinetic energy for the scalar held from the Lagrangian density is —\g tJ,v '^/ lx (j)*'W v 4 > — 
. The equation-of-state parameter w = ^ for dark energy is 


= |izhj 

V UI2 


( 5 ) 


B + vM) 

and one can see that w# < — 1 and the physically reasonable condition > 0 imply 
P(j> + P<t> = = 2 KE^ < 0, which cannot be true for a complex or real scalar held. So 

in order for < — 1 to be consistent with the positivity of —■ i-, the usual procedure is to 

hip the sign of the kinetic term in the Lagrangian density so that the dark energy density 
is — 7^2 + C(|0|) and the pressure is —^2 — V^(|0|). This is what is done in k-essence dark 
energy, and as we have mentioned, the negative kinetic term results in vacuum instability 
and difficulties in framing a viable effective held theory. We do not use this framework. 
The positivity of the kinetic energy for the usual sign of the kinetic term, 2 KE^ > 0, is 
equivalent to w$ > —1, which is the Null Energy Condition (p^ + P^ > 0) for a single perfect 
fluid with positive density^ 

In order to avoid the unwanted pathologies of phantom helds while still maintaining w$ < 
—1, perhaps the kinetic energy can be kept positive when examined under a perturbative 
departure from isotropy and homogeneity. In a universe that we know is not perfectly 
isotropic and homogeneous, perhaps describing phantom dark energy as a held theory in 
such a Oth-order framework is neither appropriate nor physically accurate. 

If we perturb the held 0(r) —> 0(r) + 50(T, r), the density and pressure are also perturbed 
(p<p(T) p^(r) + Sp^x, r), P^(r) -i> P^(r) + SP^x, r)), and we have 


_ p* + SP^ ^(|0| 2 + 0*50 + 050*) - (V + W50) 

U>eS = 


( 6 ) 


P<t> + Sp* 2^(|0| 2 + 0*50 + 050*) + {V + V'50) ’ 
where ' denotes dif ferentiation with respect to 0 and KE e ff = ^ (|0| 2 + 0*50 + 050*). In 
analogy with the Oth-order case, the positivity of the kinetic energy, 2 KE e ff > 0, is equivalent 
to ta e ff > —1, which is the Null Energy Condition (p^ + Sp$ + P^ + 5P^ > 0) for a single 
perfect fluid with positive density. 

1 Interestingly, the Null Energy Condition (NEC) is not violated for all < — 1 at present time since it 
requires Xh(Pi + Pi) > 0. If one includes both dark energy and matter components at the present time, 
utcf, can be as low as —1 — Plmo/^DEo = —1-46 and still not violate the NEC. However, dark energy would 
still have a negative kinetic term. We do not frame our discussion primarily in terms of the classical 


energy conditions of general relativity; we are merely concerned with the sign of the kinetic energy. 






However, in the perturbative approach, the equations hold separately for Oth and 1st 
orders, and the lst-order perturbation parameters are evaluated assuming the Oth-order 
equations. For a fluid model of phantom dark energy, this approach is completely consistent, 
but it is not for a phantom scalar field because of the incompatibility at the Oth order we 
saw earlier. Instead, we can define a scalar field theory for phantom dark energy in the 
background of the FLRW metric only when it is with lst-order perturbations. In that case, 
we define a field $ valid only at 1st order by 


p$(f, r) 

- P4>{t) +Sp 4> (x,r) - 2 J 2 k 2 2 +E(|$|), 

(7) 

Pf>(x,r) 

= P*(t) + SP«,(X, 

(8) 


1 (j) |2 

2KE$ = p$ + m - k 2 1$| 2 , 

(9) 


where the term proportional to k 2 is present for a field <3>(x, r) that is not spatially homo¬ 
geneous. So for an apparent value of < — 1 as measured by observational probes, it may 
be the case that w<s> = P$/p<j> > —1 and KE$ > 0, indicative of a viable scalar (real or 
complex) field theory for phantom dark energjj^J 


Cosmological Perturbation Theory 


The Friedmann equations resulting from solving Einstein’s equation for the flat FLRW 
metric (Oth order), Eq. Q, are 


n 2 

n 


8 nG 9 
^P, 



(p + 3P), 


( 10 ) 

( 11 ) 


where P = - and p and P represent the sum of the density and pressure components 
respectively. These equations lead to the evolution equation for each density component: 


p = -3 n(p + p). ( 12 ) 

2 Note that the condition for positive KE, p^ + dp^ + Pp + SP^ > 0, has terms of zeroth and first order in the 
same inequality even though first order perturbations are small compared to zeroth order terms. However, 
this inequality is especially useful in studying the sign of the kinetic energy when the first order terms 
make a non-negligible contribution to the inequality, namely when the zeroth order terms are collectively 
small. This is the case when \pp + P<p\ -C 1, or \w# + 1| <C 1. 
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Throughout this work, we consider scalar perturbations in the synchronous gauge, and 
we use the notation of [8] in the following. The perturbed metric is 

ds 2 = a 2 (r) [—dr 2 + (Sij + hij)dx l dx j ] , (13) 


and the scalar mode of h t] is written in fc-space as 


hij(x,T ) = / d 3 ke lkx \ kikjh(k,r ) + (kikj — -^-)Grj(k, r) 


(14) 


The equations resulting from solving the perturbed Einstein equation in Fourier space to 
first order are 


k 2 r]- x - nh 
k 2 f] 

h + 2'Hh — 2k 2 rj 
h + 6 rj + 2 H(h + 6 rj) — 2 k 2 r] 


AnGa 2 ST/, 

(15a) 

4TiGa 2 {p + P)6, 

(15b) 

- 8nGa 2 5Ti, 

(15c) 

— 24:7iGa 2 (p + P)a, 

(15d) 


where 9 is the divergence of the fluid velocity Uj, (p+P)c r = —(kikj — ^Sij) S* • where S® ■ is the 
anisotropic shear perturbation, and h and r) are the scalar modes of the metric perturbation. 
The stress-energy tensor is given by 


n = -(p+sp), 

A = (p + pa, 

Tj = (p + JPJiC + EC, ^, = 0. (16) 


The conservation of energy-momentum, = 0, gives (using S = dp/p) 


h 


6 = -n{l-3w)6- 


w 


1 + w 


9 + 


(g-)‘ 

(17a) 

Woet-e*. 

1 + w 

(17b) 


Eq. ( [I7| ) is valid when considering each fluid component or the total fluid, but Eq. (15) is 
valid only for the total fluid. The anisotropic shear stress is 0 (cr^ = 0) throughout, and 
in what follows, we use c — G — 1. As is typically done to define the coordinates of the 
synchronous gauge, we choose 6 C = 0, for cold dark matter. 


5 





( 18 ) 


SP/p for a given fluid component is in general given by 


— = c 2 s S + (c 2 s - c 2 a )31d(l + w)^, 
p k z 

where c s is the fluid’s sound speed and c 2 a = P/p = w + wp/p is defined as the square of the 
fluid’s adiabatic sound speed [9]. For a barotropic fluid, (? s = cl, and c 2 a = w for constant 
w. Even though dark energy can have a barotropic equation of state, treating it like an 


adiabatic fluid (for which Eq. (18) reduces to 5P = c 2 s 8p) would imply imaginary sound 
speed and instabilities in dark energy, so we use this general relation between 5P and 8p. 


Condition for Positive Kinetic Term 

We want to determine if p$ + 8p^ + P^ + dP^ > 0 (w e g > —1) using observational 
constraints. For positive density, this inequality is equivalent to 


1 + ^0 + Wcf, + SP^/pif, > 0 . 


(19) 


Using Eqs. (18) and (17b) to put 5^ in terms of 6$ and 9^, this inequality becomes 


(1+UV>) ( 1 + -jp 


where V = 9/k. 


1 + ^ dV$ 

c L da 


o, + V<p 


(l + ^)(l-3c^) 




+ 3 c S( p — ‘iWib / + 


dw^ a 
da 1 + 


(20 


For Wtf, < —1, Inequality (20) is equivalent to 


n 

~k 


1 + c^dV <p 

C U d a 


a + y<j> 


(l + ^)(l-3c^) 


S(f)J 


''scf) 


+ 3 c S( p — 3 Wd, /■ + 


dw^ a 
da 1 + 


< -1. (21) 


For the case of a scalar held [9] , 


2 _ dP* 


astf> — 


dpo 


= 1. 


rest frame 


Inequality (|2l|) then reduces to 


n 


p = 


k 


da 


+ y<p I { — 1 — 3w<p} + 


dvjff, a 
da 1 + w^J 


< - 1 . 


( 22 ) 


(23) 


Even with the constraints |V^| < 1 and \V^\ = adht d ^// 1 < 1, it is mathematically possible 
for Inequality (23) to be satisfied, but we must solve for from Eqs. (15) and © to 
determine for certain. 


> 0 , 
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Dark Energy with Constant w 


For dark energy with constant w, Eqs. (15) and (17) can be solved explicitly for different 


eras dominated by one fluid component. Assuming one fluid component with constant w 
during an era dominated by that fluid component, RL is given by 

2 


n = 


(24) 


(3w + 1 )t 

For w < —1/3, r G (—oo, 0) for a G (0, oo), whereas r is positive for w > —1/3. For working 
between different eras, we find it more convenient to work in terms of a rather than r. r 


can be related to a via Eq. (10). 


We know that the universe is homogeneous and isotropic only on large scales; on smaller 
scales, the universe obviously is not described by a FLRW metric, even one with small lst- 
order perturbations. And since we are only concerned with dark energy, which we know to 
act only on large scales at which its repulsive force is detectable relative to the attractive 
force of gravity that keeps smaller structures bound together, we do not consider large k. 
As one might expect, the assumption of lst-order perturbations (|5[ <C 1 , \SP/P\ <C 1, 
|V^| -C 1, and their derivatives with respect to r) is not always satisfied for large k. The 
spectra of Type la supernovae that have enabled the detection of universal acceleration span 
in redshift from z ~ 0.3 to z ~ 2 [10], which corresponds to spectral shifting over a range 
of 10 -4 MpW 1 < k < 8 x 10 -4 Mpc -1 . The late-time Sachs-Wolfc effect, which is evidence 
for the effect of dark energy on superclusters and supervoids in the CMB, also occurs on 
similarly large scales. An acceptable theory of dark energy must at least be valid for a 
relevant range of large scales, at least for 10 ~ 4 Mpc -1 < k < 8 x 10 -4 Mpc^ 1 , and this is 
the range of k we consider in the following analysis. 

Radiation Domination 

During radiation domination (RL = 1/r and r = a( ^4> ad with constant ^ —1 


and sound speed c S 0 is obtained from Eq. (17) and from the non-gauge mode for h for 
k << RL during the radiation era [ 8 ], h — A(fcr) 2 . We set the constants of integration for 
the decaying modes of the full solution for V™ d to 0, and this gives a real solution for V£ ad . 
For k RL, VT ad is 


t/- rad 

W 


3 3 / 2 a 3 Ac 2 S( f ) k 3 


2 n / 2 (7rp r0 ) 3 / 2 (4 + 3 c 2 S(j> - 6^) ’ 


(25) 
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where A, which is in general a function of k, is a constant of integration with respect to r 
(from the expression for h ) and p r0 = '-yyQ r0 is the present-day density of radiation, which 
we define to be at a = 1. We use H 0 = 67.3 km/s/Mpc and Q r0 = 9.24 x 10 -5 , both of 
which are consistent with Planck’s reported values [3]. 

Since the comoving curvature perturbation 1Z equals the square root of the scalar power 

/ / \ n s —1 


spectrum = ±\ A s [ y- 


for adiabatic perturbations, and since 


77 = -T 7 +^[h + 6 p] + H^ (26) 

in the synchronous gauge, we can use Planck’s constraint on A s and n s and expressions for 
r) and h to get a value for A in V£ ad . 

On distance scales for which linear perturbation theory applies (roughly k > 0.1 Mpc -1 ), 
we are agnostic concerning the sign of 1Z since it appears only in squared form in all obser¬ 
vations, including the CMB angular power spectrum and the matter power spectrum. The 
sign of A, which appears in all the expressions for super-horizon perturbations in the early 
universe, determines the sign of density and velocity perturbations, and A is determined 
by the sign of 7Z. The Press-Schecter formalism is a good approximate description of the 
formation of galaxies and clusters, and it assumes positive matter density perturbations. 
However, it operates on smaller scales for which linear perturbation theory breaks down and 
higher order terms may dominate. So we will deal with both signs in our analysis. 

From [8], for adiabatic initial conditions and k <C R during radiation domination, 


t] = 2A — 


5 + 4 R v 
6(15 + AR V ) 


A(kr ) 2 , h = A(kr) 2 , 


and 


V = (1 - R v )V^ + R v V v = --A(kr) 3 ( 1 - R v + R, 


23 + 4 R, 
15 + 4 R h 


(27) 


(28) 


where R u = - p + p is the contribution of neutrino density out of the total radiation 
density from photons and neutrinos. Using our aforementioned value of p r0 , = 

(71W/8)(4/ll) 4 / 3 , and Planck’s value of N„ = 3.046, we get R u = 0.409. Combining Eqs. 


(26), (27), and (28) gives 


A(k) = 


-FV^^(15 + 4 R v 


(29) 


2(5 + AR V ) 

We see that 1Z is constant with respect to r on super-horizon scales. We use ln(10 10 H s ) = 
3.064 ± 0.023, n s = 0.9667 ± 0.0040, and a pivot scale of fc* = 0.05 MpW 1 from Planck’s 
2015 results [3|. For most values of relevant k we consider, A TlO' 5 . 










Using Eq. (25), Ineq. (23) becomes 


3Ak 2 a 2 (c 2 .J3wa — 1) — 4) 

// =- * ——y < -l. ( 30 ) 

16np r0 ( 3c^ - + 4J 


We see that for w < —1, the inequality can only be satisfied if A > 0. 


In Ineq. (23), A is a prefactor of the entire lefthand side, so the magnitude of A largely 
controls whether the inequality is met. Without constraints on A, we can find a scenario in 


which Ineq. (30) is satisfied for some ranges of parameters, and models of dark energy with 
W 0 < —1 do not obviously have a negative kinetic term. For example, for k = 10 -4 Mpc -1 , 
= —1.1, and A = 44445, the inequality is satisfied for a > 1.9 x 10~ 4 , and \V™ d \, \V™ d \ = 
\aH V™ d '(a)|, \y d \, \8P^/P^\ < 1 throughout. However, when we constrain A for adiabatic 
initial conditions, we get that |H| is too small for Ineq. (30) to be satisfied. // is made most 
negative for biggest w$ (w^ —» —1), maximal a (a = a eq = 1/3043, where a eq is when the 
density of radiation and matter are equal), and maximal k for which k % is still accurate 
(k ~ Hiaeq) = 0.0066 Mpc^ 1 ). However, when we use these parameters with our value of 


A{k = ' H(a eg )) from Eq. (29), p ~ — 3 x 10 6 , and Ineq. (30) is clearly not satisfied. 


The assumption of isocurvature perturbations changes the contribution from perturba¬ 
tions. Isocurvature perturbations are fairly tightly constrained by observations. The maxi¬ 
mally allowed contribution from isocurvature modes a non _ a di from observational constraints 
from Planck [12] are 


cold dark matter density mode : (—1.5%, 1.9%), 
neutrino density mode : (—4.0%, 1.4%), 
neutrino velocity mode : (—2.3%, 2.4%). 


Using the maximal amount of isocurvature contribution in each case to maximize the mag¬ 
nitude of p has a very negligible effect on perturbations’ values. 

The perturbations in the radiation era for different isocurvature modes are given in several 
places in the literature. For the synchronous gauge, h, S C1 and in the cold dark matter 
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density isocurvature mode (GDI) are initialized by 

h = A C di—f==Hqt — -A C di 


S r = A, 


\ZrTpQ 
h 


f^cO^mO tt 2 2 
H 0 t 


n 


r 0 


CDI 


S') — ~7,^i 


( 31 ) 


where f2 c o and Q m o are the present-day cold dark matter and total matter density contribu¬ 
tion respectively. We can then obtain the initial condition for V£ ad for k <C Zi for the CD1 


mode using Eq. (17) with h from Eq. (31). Using a best fit from Planck’s TT spectrum and 


low -l polarization modes [12], the magnitude of the primordial power spectrum at a scale 
of k\ = 0.002 Mpffi 1 is V^x — 1-4 x 10~ n for the GDI mode and V] zn = 2.4 x 10 -9 for the 
adiabatic mode, and at k 2 = 0.100 Mpc” 1 , P| x = 4.7 x 10” 13 and = 2.1 x 10~ 9 . They 
use 

In k — In ki \ T / 2 


Vab(k) = Exp 


In k — In k 2 
In ki — In k 2 


Hv l ab ) + 


In k 2 — In k\ 


ln(P, 


ab) 


(32) 


as their parameterization for the primordial power spectrum, and a, b = TZ, Z. We can 


specify the constant Acdi using Eq. (31), Eq. (32), and the relationship between the 


entropy perturbation and the isocurvature primordial power spectrum, 

S n <U 


Sry = 


1 + W a 1 + Wr 


= Z\fV: 


xx 


(33) 


for the isocurvature mode denoted by a. After specifying Acdi from data, we can use h 
from the GDI mode to solve for V™ d in the GDI mode using Eq. (17), and then we can 
approximate for k ZL. Using the sum of the adiabatic and CD1 modes for V£ ad in fneq. 
(23), we find that the most negative /i can be is ~ —2 x 1CU 5 (for the same parameters listed 
for the purely adiabatic case above), not much less compared to the exclusively adiabatic 
contribution of ~ —3 x 10 -6 . 

For the other isocurvature modes, the initial conditions for the perturbations have the 
same functional form since they are expressed as power series, but the coefficients of their 
terms (ignoring the constant of integration specified by the isocurvature primordial power 
spectrum for each mode) of the perturbations are considerably smaller in magnitude com¬ 
pared to those of the CD1 mode, so their maximal contribution allowed by observations does 


not significantly alter the chances of satisfying Ineq. (23). 

So we conclude that phantom dark energy with positive kinetic energy is not possible 
during the radiation era for constant w with lst-order corrections. 
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Matter Domination 


During matter domination (H = 2/r and r = (~~) 1//2 )> the metric perturbation h = 


W{kr) 2 (leaving out the gauge modes) is obtained from Eqs. (15) and (17). V™ att during 
matter domination with constant w$ 7 ^ —1 and sound speed c sc j, is obtained from Eq. (17) 
and h. We set the constants of integration for the decaying modes of the full solution for 
yrnatt ^ () anc [ )-}q s gi ves a real solution for V™ att . For k <C "H, 

3 3/2 c 2 Wk 3 a 3/2 


rmatt jg 


V, 


matt 


(34) 


2 5 / 2 (p m o7r ) 5/2 (5 + 9 c 2 s</> - 15w < j ) ) ’ 

where IE is a constant of integration (from h) and p m0 = is the present-day density 

of matter. We use D m0 = 0.315, which is consistent with Planck’s reported value [3]. 


Using Eq. (23), the condition for a positive kinetic term is 

3ID(c 2 (1 + 6 w 0 ) - 5)k 2 a 


p = 


< - 1 . 


(35) 


4vrp m0 (5 + 9c 2 ^ - 15 w#) 

For < —1, we see that the inequality can only be satisfied for W < 0. Since h goes 
like to a constant with respect to time multiplied by ( kr ) 2 for adiabatic perturbations for 


both the radiation and matter eras, we see that W — A via continuity of h, and Eq. (29) 
also applies to W. For negative W, p is most negative for most negative w$ (which would 


be ~ —2 from observational constraints), largest a (which is a = ode = ( 


_ ( PmO \ 3w 


DEo3w,I, 




which is when the densities of matter and dark energy are equal and dark energy begins 
to dominate for constant w^), and largest k (k ~ 'H{(1 de))- Using these parameters along 
with W from Eq. (29), we hnd that the lefthand side of Eq. (35) is ~ 3 x 10~ 6 , so Ineq. 


is not satisfied. Although our approximation k TL breaks down near the end of 
matter domination since W^de) ~ 10 “ 4 , our conclusion still holds, and for smaller values 


of a during matter domination for which the approximation is more accurate, Ineq. (35) is 
certainly not satisfied. 

When including the maximal contribution from isocurvature (and we deal with the 
CDI mode for the reasons discussed in the previous section), we hnd W from matching 
tot (a eq , k) = V^ u {a eqi k), where V^t is the sum of the adiabatic and CDI contributions 
for V£ ad , and we use the best fit that we used in the previous section on the radiation era. 
This matching barely changes the magnitude of ID from what it was in the purely adiabatic 


case, and p changes from ~ —3 x 10 6 to ~ —1 x 10 5 , leaving Ineq. (|35l) still unsatisfied 


-5 
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Notice that the constant of integration A from radiation domination and W from matter 
domination must be oppositely signed in order for the kinetic energy to be positive during 
both periods. And we found W from the condition W = A for adiabatic perturbations, and 
this is approximately true also when isocurvature contributions, which are small, are con¬ 
sidered. So we see that positive kinetic energy during radiation domination is incompatible 
with positive kinetic energy during matter domination for k EL. 

So we conclude that positive kinetic energy is not possible during the matter era for 
constant w$ with lst-order corrections. 


Dark Energy Domination 


During dark energy domination {% = (3 ^ 2 +1)t , t = I V ? E with 
constant w$ 7 ^ —1 and sound speed c S0 is obtained from Eqs. (15a), (15b), (17a), and (17b): 


K E = S[ fpDEO 


2tt 


3-WJ. + 1 


-1 


(36) 


3^2 

where pdeo = ~^~^deo is the present-day density of dark energy. We use Ddeo = 0.685, 
which is consistent with Planck’s reported values |[3j. V E)E is given without making the 
approximation k -C EL. It has two other modes that are purely imaginary, so we set their 


constants of integration to 0 in order to have a real expression in Eq. (36). 


Inequality (23) becomes 


'2tt V ” v+1 (^£Pdeo) 1/2 , 0 , Q ^ - 3 ±*“+ ^ 
P = ~S ( — pdeo I ——^-(3 + 3 w^a 2 < 


- 1 . 


(37) 


The only way for this inequality to be satisfied is if S' < 0. S is found through matching 
ymatt( aDE ) = V EE (aoE) for 10 ~ 4 Mpc -1 < k < 8 x 10 -4 Mpc -1 and —2<w< —1, our 
allowed ranges for k and w^. We do the matching for the full solution for V™ att since k <C EL 
is not accurate at a^E- 

For either the purely adiabatic case or the mixed case with maximal isocurvature contri¬ 
bution, this matching implies 10 -10 < S < 10 -5 , and this implies — 6xlCE 4 < p < —5x ICE 8 , 


leaving Ineq. (23) unsatisfied. 


In fact, using Eq. (17b), one Ends that p = S EE . So p < — 1 implies S EE < —1, which 
would break our assumption of \SP E \ -C 1 for perturbation theory, independent of S. 
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So we conclude that positive kinetic energy for dark energy as a scalar field is not possible 
using lst-order cosmological perturbation theory for the FLRW metric and constant w<p and 


Non-Constant w$ Models 


For models with 7 ^ 0, the dynamics of fneq. (23) are different due to the term 
and w^s effect on V^. 

Dark energy with constant < — 1 has density that behaves as a power law in a as 
P(j> — Ppo a _ 3 b"* +1 ). All these models lead to a big rip [H] since p^ —* 00 at a finite time in 
the future and all bound structures are ripped apart before that time due to the increasing 
repulsive force of dark energy. We consider in this section two models with p^ that increase 
more slowly than any power law in a. We examine a model leading to a little rip |15], 16] . 
which means that p^ —> 00 as t —> 00 and all bound structures are ripped apart eventually. 
We also consider a model leading to a pseudo-rip H 7 ], which means that p^ —* as t —* 00 
and all bound structures at or below the threshold determined by p^ rip apart. 

The little rip parametrization we consider is given by 

3ct 


Plr ~ PDEO 


2 n 1/2 

^Pde 0 


In a + 1 


(38) 


This parametrization was fitted to recent supernovae data with the best fit of a = 2.26 x 
10 -6 Mpc -1 [15]. The fitting essentially ensures Wi r ~ 0 around a = 1. And note that 
pir(a = 1) = pdeo by construction. The equation-of-state parameter, obtained from Eqs. 


( 10 ) and ( 11 ), is 


Wlr = -1 


(ip'M 


= -1 


3 Pw | In a + yj pdeo/ a ' 

The pseudo-rip parametrization we consider is given by 


Ppr — Pdeo I - 


/ ln l 7+tu + 


(39) 


(40) 


AHjir + l 

where / = 10 ~ 23 and s = 48. This parametrization was also fitted to recent supernovae 
data with the best fit of b = 0.01078 H 33 . and p pr (a = 1) = pdeo by construction. The 
equation-of-state parameter is 


w P r = — 1 


a Ppr( a ) 


3 P, 


= —1 — sb- 


pr 


3(1 + o/)( 1 + ci(b + /)) In 


_J_ + 1 
A+/ + b 

a J 


(41) 
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Note that both wi r and w pr are strictly less than —1 for all a , and they approach —1 as 
a —> oo. 

We consider adiabatic perturbations during the radiation era (h = A(kr) 2 ), matter era 
(h = W(kr) 2 ), and dark energy era (h put in terms of and its derivatives using Eqs. 
(15a), ( 15b[ ), (17a), and (17b)) for the little rip and pseudo-rip parametrizations, and we 
numerically solve the differential equation we obtain for V^. When solving for a model with 
non-constant w<p, is no longer analytic in general, and one must specify initial conditions 
beyond the one parameter that is constrained by the primordial power spectrum in the 
radiation and matter eras. For reasonable initial conditions for both parametrizations (given 
in the captions of Figs. and 0 ) and A and W constrained from the primordial 

power spectrum, we plot /i in Figs. ( 0 . § , and 0 for both k = 10 4 Mpc 1 and k = 
8 x 10 -4 Mpc -1 , the upper and lower bound of observationally relevant k that we consider. 
For the radiation and matter eras, the inequality is satisfied for some of the era (where 
/i < —1), but not for the whole era; this is the case for all k, independent of the initial 
conditions chosen for VT ad and V7" att and the signs of A and W. Adding in the isocurvature 


mode for V£ ad (which uses h from Eq. (31)) and propagating this addition through the next 
eras does not change the magnitude or shape of the results in any of the three figures very 
much at all. 

In Fig. 0, we plot n from the beginning of the dark energy domination era to a = 1 
(present day). We End that one can choose a set of initial conditions such that p < —1 for all 
of k G (10 -4 Mpc -1 , 8 x 10 -4 Mpc -1 ) while still meeting the requirements of smallness for 
perturbation theory. Even though Wi r and w pr fail to provide positive kinetic energy for the 
whole of the radiation and matter eras, in theory, dark energy could be a phenomenon that is 
active only during its domination era and not before (as we have no observational evidence for 
dark energy’s existence before its era, strictly speaking). If dark energy becomes active only 
during dark energy domination as the result of, say, some spontaneous symmetry breaking, 
then there is no necessary continuity with the dark energy perturbations we calculated for 
the matter era. Fig. 0 indicates that it may be possible to have phantom dark energy 
that satisfies a held theory with positive kinetic energy for all relevant observational k and 
a (that is, up to the present). However, in the future (a > 1) for the models we examine, 
H generally becomes > —1. Also, one can be more precise by considering an era dominated 
by both matter and dark energy. 
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FIG. 1. We plot p, the lefthand side of Ineq. (23), for k = 10 -4 Mpc~ 4 (blue lines) and k = 8 x 
10” 4 Mpc -1 (red lines) with chosen initial conditions VJ ad (10~ 5 ) = 10~ 2 and VJ arf (10 -5 ) = 2 x 10 -4 
during the radiation era for adiabatic perturbations for both the little rip (solid lines) and pseudo¬ 


rip (dashed lines) parametrizations from Eqs. (39) and (41). All perturbations are sufficiently 
small. The top set of patterned lines is red. (Color plots available in the online version of this 
work.) We see that the condition for positive kinetic energy can be satisfied for some but not all 
of the radiation era. 


DM-DE Domination Era 


We also consider an era with both dark matter and dark energy, EL = cl(^-(pde{o) + 
Pco gC 3 )) 1 / 2 , where Pde{o) is the density of dark energy for a generic w^a) (constant or 
not), and p c o = '-^-D c q is the present-day density of dark matter. We use n c0 = 0.266, 
which is consistent with Planck’s reported values [3]. Because 6 C = 0 in the synchronous 


gauge and 5P C = 0, Eqs. (15) and (17) are identical to the ones used for the dark energy 


era except for Eq. (15a), which contains contributions from both S c and 8 ^ on the righthand 


side. We use Eqs. (15a), (15b), (17a), and (17b) to obtain our differential equation in terms 


of V<f) and its derivatives. 

We have already seen that positivity of kinetic energy is not satisfied for constant w$ 
during the matter and dark energy eras, so one would expect a similar outcome for the 


DM-DE era. Indeed, one can match V™ d (a eq ) = V EM DE (a eq ) along with the derivatives to 
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FIG. 2. We plot /x, the lefthand side of Ineq. (23), for k = 10 4 Mpc 1 (blue lines) and k = 


-3 


8 x 10 4 Mpc 1 (red lines) with chosen initial conditions V™ att (a eq ) = 0.4 and V™ att (a eq ) = 2x 10 
during the matter era for adiabatic perturbations for both the little rip (solid lines) and pseudo-rip 


(dashed lines) parametrizations from Eqs. (39) and (41). All perturbations are sufficiently small. 
The top set of patterned lines is red. (Color plots available in the online version of this work.) We 
see that the condition for positive kinetic energy can be satisfied for some but not all of the matter 


era. 


set the initial conditions for solving the differential equation, and we find the /i is not < — 1 
by a wide margin. 


However, if dark energy is considered active only from a^E onward, then y EE ~ DM j s not 
necessarily constrained from continuity with and its derivatives before this era, leaving 
free the choice of initial conditions. Along with trying to achieve /x < — 1 for a G (a^B, 1) 
and k G (10 -4 Mpc -1 , 8 x 10 -4 Mpc -1 ), one has to satisfy the smallness of dark energy 
perturbations and the smallness of S c , the latter of which can be checked via the use of Eq. 
(15a). After scanning the parameter space of initial conditions for y E ~ DE ) we found that 
the positivity of kinetic energy is not satisfied for the relevant ranges of both a and k , even 
if it can be satisfied for some a and k. And even for initial conditions that cause /x < — 1 for 
all a G (acre, 1) for a given £;, /x still grows to be > —1 for a > 1, and |5 C | is often > 1 for 
other values of k with the same initial conditions. In Fig. Q, we plot /x for W[ r , w pri and 
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FIG. 3. We plot n, the lefthand side of Ineq. (23), for k = 10 -4 Mpc -1 (blue lines) and k = 
8x 10 -4 Mpc -1 (red lines). We chose V EE (clde) = 0.6, V EE (clde) = —3x 10 -4 , and V EE (clde) = 
—8x 10 -8 during the dark energy era for the little rip parametrization (solid lines) and V EE (clde) = 
0.6, V EE (clde) = — 6 x ICG 4 , and V EE (clde) = —2 x 10“' for the pseudo-rip parametrization 
(dashed lines). All perturbations are sufficiently small. The top set of patterned lines is red. 
(Color plots available in the online version of this work.) We chose initial conditions such that 
the little rip parametrization does not have fi < — 1 for all k while the pseudo-rip parametrization 
does. We see that the condition for positive kinetic energy can (apparently) be satisfied for all of 
the dark energy era, but see the text for further details. 


constant for sample initial conditions. 


Dark Energy with > — 1 


It is interesting to note that for the case of > —1, the only change in Ineqs. (21) and 
(23) is < —> >, and the magnitudes of fi of our best attempts to satisfy positivity of 


kinetic energy for constant w$ does not increase substantially with a slight change in the 
magnitude of w# from a value less than —1 to one bigger than — 1 . Since V<f, was computed 
in different eras for a generic w </, ^ —1, they apply for both < — 1 and w$ > —1. 
So for > —1, it is easy to see that /i > — 1 for either choice of sign for the constants of 
integration during all single-component eras, meaning that negative kinetic energy is avoided 
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FIG. 4. We plot n, the lefthand side of Ineq. (23), for k = 10“ 4 Mpc" 1 (blue lines) and k = 
8 x 10~ 4 Mpc" 1 (red lines) during the DM-DE era for the little rip parametrization (solid lines), 
pseudo-rip parametrization (dashed lines), and constant w$ = —1.1 (dot-dashed lines) for a G 
(0.61,1), where 0.61 is close to clde for all 3 parametrizations. We chose the following initial 
conditions for k = 10" 4 Mpc" 1 : little rip: V EM ~ DE (0.61) = 0.24, V em ~ de (0.61) = -1.0 x 10~ 4 , 
yDM-DE ( 0-61 ) = _ 2 .9 x 10" 4 , i/'J M " DS (0.61) = 5.3 X 10" 3 ; pseudo-rip: V EM ~ DE (0.61) = 0.22, 

; (0.61) = 1.8x 10“ 4 , vf 


yDM-DE = —l.lx10~ 4 , V em ~ de (0M) = 1.8x 10“ 4 , V% M ~ DE {0.61) = 1.3x 10~ 3 ; constant 


w 4>- V, 


DM-DE 


(0.61) = -7.3 x 10" 2 , V em ~ de {0.61) = -4.0 x 10" 5 , V em ~ de {0M) = 4.3 x 10" 5 , 


V EM ^(0.61) = 9.7 x 10 4 . We choose the following initial conditions for k = 8 x 10 4 Mpc 4 : 
little rip: V EM ~ DE (0.61 ) = 0.0, V EM ~ DE {0 .61) = -4.5 x 10" 4 , V em ~ de {0M) = 8.0 x 10" 4 , 
yDM-DE ^ 0 61 ) = 2 ? x 10 -2. p Seu do-rip: V EM ~ DE (0.6l) = 0.0, V EM ~ DE {0.61) = -5.3 x 10” 4 , 
yDM-DE ( 0 61 ) = 2.7x 10" 3 , V° M ~ DE (0.61) = 1.4x 10" 2 ; constant W(j> : V EM ~ DE (0.61 ) = -7.2 x 
10 - 2 , yDM-DE^ o 61 ) = —3.6 x 10" 5 , V em ~ de (0M) = 1.2 x 10" 4 , V EM ~ DE (0.61) = 9.4 x 10" 4 . 
All perturbations are sufficiently small. The top set of patterned lines is red. (Color plots available 
in the online version of this work.) We see that the condition for positive kinetic energy can be 
satisfied for some but not all of the DM-DE era. 


for constant w$ > —1. However, for constant w$ > —1 in the DM-DE era, we have checked 
that n can be < — 1 (negative kinetic energy for > — 1) for certain scales and times with 
the right choice of initial conditions while the smallness of perturbation theory is satisfied, 
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FIG. 5. We plot /i, the lefthand side of Ineq. (23), for k = 10~ 4 Mpc -1 (blue line) and k = 
8 x 1CT 4 Mpc -1 (red line) during the DM-DE era for constant = —0.99 (dot-dashed lines) 
for a £ (0.61,1) and for the same initial conditions for the constant w$ cases in Fig. (|4j). All 
perturbations are sufficiently small. The top patterned line is red. (Color plots available in the 
online version of this work.) We see that positive kinetic energy is not satisfied for all of the DM-DE 
era. 


just as we analogously found /i < — 1 (positive kinetic energy for w^ < —1) possible in 
this era for constant < —1 for certain scales and times! See Fig. (J5]) for an example 
of negative kinetic energy from a > — 1 model. Note that the constant w^ lines from 
Figs. Q and (J5]) are very similar. We suspect that non-constant parametrizations for non¬ 
phantom dark energy can also give rise to particular scales and times during which kinetic 
energy is negative. This prospect of quintessence dark energy’s violation of the positivity of 
kinetic energy is intriguing and should be pursued in future research^] 


Conclusion 

In this work, we have studied models of phantom dark energy completely within the 
confines of general relativity to discern whether or not they can be represented as a scalar 
3 For the case of = —1, SP^/= c^S^, and Ineq. (19) becomes 1 + ^(1 + c^) +w ;0 >0. If we multiply 


every term in Eq. (17b) by (1 + w), we see that 6^ = 0 for = 0, and Ineq. (19) becomes l + u >0 = 0 > 0. 
So negative kinetic energy is avoided for wyj = — 1 . 
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field theory with positive kinetic energy. While it is clear for FLRW space at Oth order that 
every phantom model of dark energy fails to give a positive kinetic term for a scalar held 
theory, it is not as clear at 1st order. For a scalar held defined for phantom dark energy 
for the more observationally accurate lst-order FLRW space and not for Oth order, it is 
mathematically possible to satisfy positivity of kinetic energy (// < — 1) for certain relevant 
ranges of a and k. 

For constant w<f> models during the radiation and matter eras constrained by the primor¬ 
dial power spectrum, positivity of kinetic energy is not satished, and for the dark energy 
era, perturbation theory would have to be violated in order for fi < — 1 since /i = 8 EE . For 
the DM-DE era, p > — 1 if y EM ~ DE j s matched up with the observationally constrained 
VT ad for models with constant w^, and positivity of kinetic energy is not satished. 

For the models with non-constant w$ we consider during the radiation and matter eras, 
despite the freedom to choose initial conditions while still meeting constraints from the 
primordial power spectrum, we find that positivity of kinetic energy is satished only for 
limited ranges of a and k. During the dark energy era, our models with non-constant w$ 
can satisfy positivity of kinetic energy for all relevant a and observational k given the right 
initial conditions. But generally fi > — 1 for the future (a > 1), and the more accurate 
description of the DM-DE era provides different results. For our models with constant or 
non-constant during the DM-DE era, positivity of kinetic energy can be satished only for 
some a and k without violating perturbation theory’s assumption of smallness, and /i > — 1 


in the future. In fact, one could be even more accurate and assume contributions to Eq. (15) 
and H from radiation, matter, and dark energy, and the added constraints on the dark energy 
perturbations from the requirement of smallness for radiation and matter perturbations may 
make the satisfaction of the positivity of kinetic energy even more difficult. 

For the case of > —1, we saw that violation of the positivity of kinetic energy is 
possible and should be further examined. We suspect that the more accurate assumption of 


contributions to Eqs. (15) and from radiation, matter, and dark energy and constraints 
from the smallness of all the components’ perturbations would make the satisfaction of 
negative kinetic energy for non-phantom parametrizations more difficult. 

So overall, in the context of general relativity, we fold that phantom dark energy does 
not provide a scalar held theory with positive kinetic energy in lst-order FLRW space 
for all relevant scales and times, at least for const ant-iu models and the others we have 
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tested, but we do point out that phantom dark energy is not categorically indicative of 
negative kinetic energy for a scalar held theory, and, as far as we know, this work is the 
first attempt in explicitly checking the possibility of saving phantom dark energy as a scalar 
held theory purely through perturbations. By the same token, we found that quintessence 
models may have negative kinetic energy for certain times and scales, and this should be 
further explored in future work. Other scenarios in the framework of FLRW cosmology 
can be studied rigorously for the positivity of the kinetic energy of dark energy which may 
serve to more precisely quantify the effects we demonstrate. Higher-order perturbations and 
the back-reaction of perturbations on the background would be an important consideration. 
Additionally, if dark energy is not a perfect fluid, different conditions may apply for the 
positivity of its kinetic energy. And if dark energy and dark matter are coupled, the condition 
for positive kinetic energy changes, although some models may suffer from instabilities during 
the early universe, as |9] points out. For the covariant model they propose, perhaps the 
divergence in the gauge-invariant curvature perturbation can be pushed back early enough 
to a time for which quantum gravity dominates the relevant length scales so that there is no 
divergence for the classical regime. Other frameworks (such as modified gravity, the presence 
of quantum effects, vector field dark energy, non-perturbative effects) and metrics (such as 
Bianchi and Tolman-Bondi models) for a universe deviating from a perfectly isotropic and 
homogeneous one can be studied to determine the viability of phantom dark energy. Along 
with these different frameworks come different manifestations of observational constraints. 
We leave these other avenues to future work. 
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